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The effects of spin-1 impurities introduced in dimerized antiferromagnetic spin chains are studied 
using Lanczos diagonalization and quantum Monte Carlo simulations. Firstly, the magnitude of 
the interaction between a Ni-ion replacing a Cu-ion is derived for the case study system Cu-O-Cu. 
Then, it is shown that the introduction of S = 1 impurities leads to the suppression of the original 
spin gap of the dimerized system and to the enhancement of the antiferromagnetic correlations 
close to the impurity. These features are interpreted in terms of a valence bond picture of the 
ground state of both the pure and the doped system. The similitudes and differences between the 
underlying mechanisms that produce these effects in the cases of nonmagnetic and S = 1 impurities 
are analyzed. It is stressed that a singlet formed between the spin-1 impurity site and the two spin- 
1/2 neighbors spins behaves as an effective nonmagnetic vacancy which may be partially responsible 
for the enhancement of the antiferromagnetic correlations on the whole system. Implications of 
these results for experiment are also discussed. 
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I. INTRODUCTION 

One-dimensional or quasi-one-dimensional magnetic 
systems have several fascinating properties which con- 
tinue to attract an intense theoretical activity. One of 
these properties is the presence of a spin gap in anti- 
ferromagnetic Heisenberg chains with integer spin [QfJ] 
and with spin-1/2 in ladder geometries |i^J], as well as in 
chains with frustrating interactions ||, and in the pres- 
ence of dimerization. In the latter systems, the dimer- 
ization frequently appears as a consequence of the spin- 
phonon coupling, i.e. it is due to a spin-Peierls (SP) in- 
stability. The interest in the spin-Peierls phenomena was 
recently revived after the first inorganic SP compound, 
CuGe03, was found (6). A considerable theoretical ef- 
fort has been devoted to understanding the properties of 
this compound @,||. Most of these studies have been 
based on a one-dimensional (ID) dimerized Heisenberg 
model, i.e. retaining only the spin degrees of freedom and 
neglecting the spin-phonon dynamics. Although most 
dimerized systems are spin-Peierls systems, there are also 
several compounds where the intrinsic structure leads to 
dimerized spin chains. Actually Cu(N03)2 ■ 2.5H2O was 
the first inorganic alternating-chain system to be inves- 
tigated More recently, inelastic neutron scattering 
measurements and susceptibility data indicate that the 
compound C11WO4 consists of weakly coupled alternat- 
ing chains |10|. Another example of structural dimeriza- 



tion is provided by the Cu02 chains which are present 
in Sri4Cu2404i (TlJ which also contains CU2O3 ladders. 
In addition, the compound (VO)2P207, which was previ- 
ously believed to be composed by two-leg spin ladders, is 
actually better described by alternating spin chains [12| . 
The experimental introduction of impurities in these 



magnetic systems has provided additional information in 
the quest to examine the nature of their ground state 
and excitations In particular, inorganic systems, 

such as CuGe03, can be easily doped with both non- 
magnetic and magnetic impurities i.e. by partially re- 
placing the Cu 2+ ions by Zn 2+ and Ni 2+ ions, respec- 
tively. The effect of impurities on CuGe03 has been 
investigated using neutron scattering [ fj^]l5| , magnetic 
susceptibility measurements Q,^6|, specific heat mea- 
surements ]l7| , and spin electron resonance [Q. The 
main results of these studies are the collapse of the spin 
gap Jl4||l6| , and the competition between the SP and the 
Neel states 14 17j], when a small amount of impurities 
is introduced. From the theoretical point of view, the 
introduction of non-magnetic impurities, which for the 
nearest-neighbors Heisenberg model leads to the appear- 
ance of chains with open boundary conditions (OBC), has 
been analyzed using conformal field theory and Monte 
Carlo simulations jl9],^0| . The main result of these stud- 
ies is the appearance of enhanced staggered spin order, 
which is maximized near the vacancies. More recently, 
it has been shown that spin-1/2 states are localized near 
vacancies in dimerized chains and other systems [^l]-^3) . 
The experimentally observed stabilization of a Neel state 
upon Zn doping in dimerized chains can be related to 
the enhancement of antiferromagnetic (AF) correlations 
near the vacancies which arises due to the presence of 
these bound states. The effect has been studied theo- 
retically using both analytical |24| and numerical tech- 
niques pl] , p2]| . These bound states are expected to pro- 
duce spectral weight inside the original gap p5| j, and 
there is already direct experimental evidence compatible 
with the existence of these low-energy excitations [^6],^7j . 
Results very similar to those obtained in dimerized chains 
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have also been obtained in ladder systems both experi- 
mentally H| and theoretically 

In the present article, the previous studies per- 
formed on ID AF spin-1/2 dimerized Heisenberg 
chains J2l],^,^5| are extended to the case of magnetic 
impurities, in particular spin-1 impurities. In previ- 
ous literature, it has been studied the related case of 
spin-1 AF Heisenberg chains doped with non-magnetic 
impurities which become, again, equivalent to chains 
with OBC |]30| , ^l| . The most important prediction in 
this context is the presence of isolated spin-1/2 states 
on the chain edges. Experimentally, magnetization 
and electron spin resonance studies on the compound 
Ni(C3HioN 2 )2N 3 (C104) have confirmed the existence of 
these end-chain spin-1/2 spins J52|. The end-states are 
easily explained by the valence bond solid (VBS) picture 
of the ground state of spin-1 Haldane chains |33|]. In the 
same spirit, an unifying picture based on a short-range 
resonant-valence-bond (RVB) scenario |34||35| has been 
formulated to explain the existence of spin-1/2 bound 
states in gapped spin-1/2 systems, as discussed in the 
previous paragraph pl] , p2fl . 

In the present study of spin-1/2 dimerized chains 
doped with spin-1 impurities, both RVB and VBS pic- 
tures will be used to interpret our results which were 
obtained using Lanczos diagonalization and quantum 
Monte Carlo (QMC) simulations. Our main conclusion 
is that the overall effects, such as the destruction of the 
spin gap by low-energy excitations due to bound states 
near the impurities, are similar to those generated by 
nonmagnetic impurities in spite of the fact that the mi- 
croscopic origin and details of these bound states are 
somewhat different. This difference manifests in a more 
localized enhancement of the AF correlations in the im- 
mediate vicinity of the 5=1 impurity than in the case 
of nonmagnetic impurities. However, at distances be- 
yond one lattice spacing away from the 5=1 impurity 
the spin-spin correlations behave very similarly to the 
case of doped nonmagnetic impurities. Thus, globally 
magnetic and nonmagnetic impurities apparently lead to 
similar physics including enhanced antiferromagnetic cor- 
relations. 

We wish to emphasize that in this article we are mainly 
concerned with dimerized chains, containing only spin 
degrees of freedom, and which can be realized in nature 
independently of the spin-Peierls mechanism. The spin- 
Peierls transitions recently found in inorganic materials 
only add an extra motivation to study dimerized chains 
since they can describe reasonably well certain features 
of these materials at very low temperatures. To properly 
describe some other features, the inclusion of phononic 
degrees of freedom is strictly necessary at least at the 
mean field level p^-|38[| (see e.g., discussion in Section 
IV). In this sense, our purpose is to give just indications 
that some of the results we found could be related to 
some experimental findings in spin-Peierls compounds. It 
is quite likely that similar experimental studies are going 
to be carried on in the near future on other materials 



containing dimerized chains. 

The paper is organized as follows. In Section || the 
strength of the effective Heisenberg interaction between 
a Cu-ion and an impurity Ni-ion is computed sta rting 
from a multiband Hubbard model. In Section III, the 
dynamic spin structure factor is studied in order to ana- 
lyze the suppression of the spin gap in a dimerized chain 
by the introduction of spin-1 impurities. The main differ- 
ences with respect to the case of nonmagnetic impurities 
are studied. The spin-spin correlations near the impuri- 
ties have also been investigated in order to evaluate the 
strength of the AF order enhancement. In Section IV, re- 



sults for the magnetic susceptibility obtained by quantum 
Monte Carlo simulations are presented. Finally, in the 
Conclusions section, the most important results are dis- 
cussed, together with their interpretation within a simple 
valence bond picture. Experimental consequences of our 
calculations are here also discussed. 



II. EFFECTIVE INTERACTION 

Before addressing the spin chain model and computa- 
tional techniques used in this paper, it is important to 
estimate the coupling between the 5 = 1 impurity and 
neighboring 5=1/2 ions. In previous studies based on 
Zn-doping, this calculation was not necessary since Zn 2+ 
has zero spin. However, for the case of Ni 2+ impurities 
with 5 = 1, a Heisenberg coupling between Ni and Cu 
should be part of any model that attempts to describe 
their behavior. In this section a rough estimation of the 
nickel-copper exchange will be made. The calculation 
will be carried out for the case of a 180° Cu-O-Ni bond, 
which is relevant for a variety of materials such as the ID 
compounds Sr 2 Cu03 and SrCuC>2, the Cu-based ladder 
systems mentioned in the previous section, the High-Tc 
cuprates, and several others. The results are intended to 
provide just qualitative information about the strength of 
the exchange. In principle, they do not apply directly to 
systems such as CuGe03 where the Cu-ions are linked by 
two oxygens with ~ 90° bonds. Carrying out a calcula- 
tion for this particular material would be difficult because 
small exchange couplings are expected to emerge from a 
problem involving large energy scales of the order of eV, 
and the rough one-bond calculation discussed below may 
not be enough to produce reliable results. 

The Hamiltonian analyzed in this section is given by 
H = H t + Hy + Hh , with the hopping term defined as 

(7 

- y^Jmoa{c' Nia ^co a , a + h.c), (1) 

where a denotes the spin, and a = 1,2 the two active 
orbitals of Ni (labeled as Nil and Ni2). cc u ,<j and CMi a ,<y 
are destruction operators for fermions on the Cu and Nia 
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orbitals, respectively. tcuO {tNiOa) is the hopping am- 
plitude for fermions to move between Cu and (Nia and 
O). The rest of the notation is standard. The interaction 
Hy is defined as 

Hv = Uc u ncu,in>Cu,l + U no,\n ,i 

a 

— A-c u ncu — Ajvi 22 nNia > (2) 

OL 

where n,cu,a is the number operator at the Cu site with 
spin a, ncu = ncu,1 JrTl Cu,{, and the notation for Ni and 
O follows analogous definitions. Uc u ,Uo,UNi regulate 
the on-site Coulomb repulsion for Cu, O, and Ni, respec- 
tively. Acu is the energy needed to move a fermion from 
Cu to O, while A^i is a similar quantity but to move 
fermions from the two Ni orbitals (assumed degenerate) 
to O. Finally, 

H H = - JffSi • S 2 (3) 

is the ferromagnetic Hund coupling (J# > 0) between 
the two active orbitals at the Ni site (with S Q being the 
spin at orbital a). 

This Hamiltonian can be solved exactly for the case 
of interest that contains 3 electrons, two of which are 
expected to be mostly at the Ni-ion and the other at 
the Cu-ion. From the difference in energy between the 
ground state and the first excited state with spin f in the 
spectrum the exchange coupling can be calculated. The 
matrix that needs to be diagonalized is 24 x 24. 

The parameters corresponding to the Cu and O ions 
are approximately known from previous theoretical and 
experimental work in the context of high-temperature su- 
perconductors. The actual parameter values used here 
are f§ A Cu ~ 3.6eV, t Cu o ~ 1.3eV, U ~ 4.6eV, and 
Ucu ~ 10.5eV. The rest of the parameters related with 
Ni-oxides were discussed in previous work by two of the 
present authors and collaborators f40]| , in the context of 
hole doped 5 = 1 chains. In that publication it was ar- 
gued that A^i ~ 6.06^, Um ~ 9.5eV, tmoi ~ 1.3eV, 
tNi02 ~ 0.75eV, and Jh ~ 1.3eV, based on a variety of 
experimental work and some cluster calculations. These 
will be the parameters used in the present study. 

As a warm-up, the similar cases of 180° Cu-O-Cu and 
Ni-O-Ni bonds were studied using the same method, since 
the exchanges for Cu02 and Ni02 planes are known ex- 
perimentally and, thus, a comparison between the com- 
putational results and properties of the real materials 
can be carried out. The Hamiltonian parameters used 
are the same as given before. For the Cu-O-Cu (Ni-O- 
Ni) bond the size of the matrix to diagonalize is 9 x 9 
(100 x 100) for the case of two (four) fermions in the 
problem. For Ni-O-Ni, four fermions are needed to repre- 
sent the two expected spin-1 states in this bond. For the 
copper oxides, the exact diagonalization of the Hamilto- 
nian matrix obtained from Eqs.(l-3) leads to an exchange 



Jcu-o-Cu — (100±30)meV, while for the nickel oxides it 
is Jni-o-Ni = (20±10)rneV^. The error bars are crudely 
estimated as follows: 1. the many parameters entering 
the Hamiltonian were arbitrarily changed one by one by 
±20% (to simulate uncertainties in their experimental 
determination), and in each case the new value of J was 
calculated, 2. a histogram with the distribution of J's 
was built, and from a fit of these results with a Gaussian 
function, the error bars were determined as the width of 
such a Gaussian. It was found that the parameters that 
affect the exchange the most are the hopping amplitudes 
and Acu- The Coulomb energies and Hund coupling are 
large enough such that small variations in their values do 
not affect quantitatively the results. 

Considering the roughness of the calculation, it is grat- 
ifying that the theoretically calculated Heisenberg cou- 
plings presented above are in very good agreement with 
the experimental results for both materials (see for in- 
stance Ref. 0)). Thus, it is expected that the theo- 
retical prediction for Jni-o-Cu will be at least qualita- 
tively correct as well. The result obtained in this paper 
is Jni-o-Cu — (100 ± 50)meV i.e. very similar to the 
case of copper-oxides. Then, it is here estimated that 
using a spin model with the same exchange J among the 
5 = 1/2 spins, and among these spins and the 5=1 im- 
purity, represents a good approximation to the problem 
investigated in this paper. 

III. DYNAMICAL STRUCTURE FACTOR AND 
SPIN-SPIN CORRELATIONS. 

The one-dimensional model which describes a spin- 1/2 
Heisenberg chain in the presence of spin-1 impurities is 
defined as: 

n = j 0- + 6 ^ ■ s *+! 

io 

+ {l + S io )s io -S io+ i}, (4) 

where S, denote the spin-1/2 degrees of freedom, and 
Si are the spin-1 operators located at the impurity sites 
io. The dimerization is staggered and it corresponds 
to Si = (— iyS. Periodic boundary conditions are im- 
posed unless otherwise stated. According to the re- 
sults of the previous section, in the following the case 
J' = J will be considered. Note that in Hamilto- 
nian Eq.Q) an interaction when two spin-1 impurities 
are located next to each other has not been included 
due to the result of the previous section suggesting that 

Jm-O-Ni "C Jcu-O-Cu, Jcu-O-Ni- 

Using Lanczos diagonalization and the standard con- 
tinued fraction formalism f39f| , the zero temperature 
spin-spin correlations and dynamical structure factor for 
model Eq.(4) have been computed. In order to analyze 
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finite-temperature properties a world-line Monte Carlo 
algorithm |Q suitable for the present problem has been 
implemented. This method has allowed us to reach low 
enough temperatures such that the ground state behavior 
of spin-spin correlations can be calculated. Thus, with 
this technique it is possible to treat longer chains than 
those accessible to the Lanczos method, having only a 
small finite-temperature contamination. 



A. One impurity 

Let us first investigate whether a suppression of the 
spin gap occurs when S—l impurities are introduced in a 
dimerized chain, similarly as it occurs in Zn-doped sys- 
tems. The anal ysis here follows closely the study per- 
formed in Ref. p5| for the case of nonmagnetic impu- 
rities. Since (i) the spin gap opens at k = n, and (ii) 
the most intense peak of the dynamical structure func- 
tion S(k,u) also appears at k = it, in agreement with 
neutron scattering experiments, then the present analy- 
sis starts by studying the behavior of S(tt,lj) upon the 
introduction of S — 1 impurities. 

The dynamical structure factor at zero temperature is 
defined by: 



S z 



(5) 



where S z {k) = (l/N) £ m S^e ikm , N is the number of 
sites, m labels the sites of the lattice, and the rest of the 
notation is standard. In the case of total dimerization, 
i.e. 5 = 1, it can be easily deduced that the introduction 
of a single impurity generates a peak in S(tt, to) inside the 
gap, actually at u> = 0. This is due to the fact that the 
matrix element M = (r] \S z (k)\r] ) = (1/3) (2 - (l/2)e ik ) 
is nonzero. Here \rjo) is the ground state of the 2-site 
dimer that carries the impurity i.e. \r]o) = ^7j(|0)T) — 

•\/2|l, J,)), where 1, 0, —1 label the states of the spin-1, and 
|, J, correspond to the spin-1/2. This simple mechanism is 
also present for impurities of higher spin such as Co, Jig] 
and actually includes Zn as a special case. The important 
point here is that in the bond being doped a state with a 
total nonzero spin appears after a spin-1/2 is replaced by 
a spin different from 1/2. This nonzero spin in the dimer 
where the doping occurs is decoupled from the rest of the 
lattice in the extreme limit 5—1 and, thus, produces 
weight at zero energy in the dynamical spin structure 
factor. 



B. Two Impurities 

Let us now consider the effect of two impurities located 
at arbitrary distances along a chain. If two impurities are 
located in a periodic chain of N sites, it is also easy to 
understand in the case of total dimerization 6 — 1 (and 



remembering that the ground state is selected to have a 
total spin projection on the z-axis equal to zero to mimic 
the case of 5 < 1) that if the first impurity is located at 
site 1 then a spectral weight at u — appears only if 
the second impurity is located at sites I2 which are even. 
For instance, if I2 is 2 then two free spin-1 are located in 
the same dimer which naturally will produce weight at 
zero frequency (remember that the Ni-Ni exchange has 
been neglected here since it is much smaller than the 
Cu-Cu or Cu-Ni exchanges). However, if l 2 = 3, then the 
matrix element M mentioned in the previous paragraph 
will have opposite signs for the dimers involving sites 1-2 
and 3-4 (note again that the total spin of the global 
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FIG. 1. Dynamical structure factor at k = n obtained by 
exact diagonalization using a 18 site chain with two 5=1 
impurities, total S z = 0, as a function of the position of the 
second impurity and for several values of 5. One of the im- 
purities is located at site 1, and the position of the second 
impurity is indicated increasing from top to bottom. 



state is zero), and thus the weight of the zero-energy peak 
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cancels. A similar argument holds for I2 equal to any 
odd-integer. On the other hand, for I2 even-integers the 
matrix elements of the two doped bonds contribute with 
the same sign and a low-energy peak should appear in 
the spectrum. These results deduced at large 5 actually 
are qualitatively similar when S is reduced from 1.0, as it 
can be observed in Fig. [j] where S zz (n 1 u>) obtained for 
N = 18 is shown. In this figure, it is clear that even for 
S as small as 0.05 a substantial weight inside the gap is 
present only when the second impurity is located in even- 
sites, in agreement with the discussion of the previous 
paragraph. In this figure, the bond between the sites 1 
and 2 is "strong" (Jj = J(l + 5)). As S decreases the 
positions of the impurity peaks inside the gap shift to 
larger values of lu and increase their weight, as expected. 

In Fig. § S zz {n,u) for N — 18 without impurities 
and in the presence of two 5 = 1 impurities for several 
values of S is presented. For this result an average over all 
the impurity positions has been carried out. In this figure 
the corresponding quantity for nonmagnetic impurities is 
also included for comparison. It can be observed that the 
weight which appears inside the original gap of the pure 
system is consistingly larger for 5 = 1 impurities than 
for nonmagnetic impurities. This weight inside the gap in 
the case of 5 = 1 impurities may correspond to the bound 
states formed between the spin-1, considered as two 5 = 
1/2 spins as in the VBS picture, and the two neighboring 
spin- 1/2 spins. It is also reasonable to assume that the 
weight of these bound states in the ground state is larger 
than the weight of the bound states due to the presence 
of OBC, at least in the RVB picture of the ground state 
discussed below, thus explaining the behavior shown in 
Fig. g 
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FIG. 2. Dynamical structure factor at k — tt obtained by 
exact diagonalization in the 18 site chain for the pure system, 
with two 5 = 1 impurities, and with two nonmagnetic impu- 
rities, total S z = 0, for several values of 8. These curves have 
been obtained by averaging over all impurity positions. 



C. Enhanced Spin Correlations Near Ni 

It is also instructive to analyze from a microscopic 
point of view the effects of a spin-1 impurity on the lo- 
cal magnetic order. This can be achieved by studying 
the spin-spin correlations around the impurity defined as 
S(l) =< S z (ii)S z (i\ + I) >, where i\ is the site next to 
the impurity located at i$, and I <C N. A normaliza- 
tion such that 5(0) = 1 was adopted. To predict the 
behavior of this correlation it is possible to invoke the 
short-range RVB picture of the ground state of spin-1/2 
ID chains previously used to describe the effect of non- 
magnetic impurities [2^j25|| . In this picture, the ground 
state is assumed to be dominated by the superposition 
of the two possible dimerized states, one with singlets on 
the odd bonds and the other with singlets on the even 
bonds. In the case of nonmagnetic impurities, a sim- 
ple argument shows that the correlation S(l) is enhanced 
near the impurities with respect to the same correlations 
for the pure system. The idea is simply that without va- 
cancies a given spin spends roughly half the time coupled 
into singlets with each of its two neighbors, while when 
the spin being studied is located next to a vacancy only 
one neighbor is available for singlet formation and natu- 
rally the correlation between these two spins is enhanced. 
The argument can be extended to distances larger than 
one lattice spacing easily plj . 

By using this same RVB picture for the case of spin- 
1 impurities one obtains, however, a somewhat different 
behavior. In the first place let us observe that the en- 
ergy of a dimer formed by a spin-1 and a spin-1/2 is — J,-, 
whereas the energy of a "pure" dimer made out of two 
spin-1/2 has energy — 3Ji/4. Then, the spin-1/2 next to 
the impurity should have a stronger tendency to form a 
singlet with the impurity rather than with other spin- 
1/2 sites, simply to minimize the energy. In addition, 
and most importantly, is that the ground state energy 
of the three site cluster formed by the impurity and its 
two spin-1/2 neighbors is equal to (— 2J), independent of 
5. If one considers a subsystem with four sites contain- 
ing the 5=1 impurity, it can be deduced that in order 
to minimize the energy this subsystem prefers to form a 
three-site cluster leaving a spin-1/2 unpaired rather than 
two dimers, one of them containing the impurity. This is 
true at least for 5 < 1/7 if the two dimers are formed on 
"strong" bonds, and for all 5 > if they correspond to 
"weak" bonds. Thus, it is expected that states contain- 
ing the ground state of this three-site 1/2-1-1/2 system, 
which corresponds to a singlet, plus an unpaired spin 
and a RVB configuration for the rest of the chain will 
have a large overlap with the actual ground state of the 
problem. The coefficients of the states in the three-site 
system ground state favors the AF correlations between 
the three sites. These two factors imply that: 1. there 
is a strong AF correlation between a site next to the im- 
purity and the impurity itself, 5j mp (l), 2. there is an 
enhancement of the correlation between both sites next 
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to the impurity, S' cross (2), and 3. there is a reduction 
of the correlations between site i\ (the site next to the 
impurity) and the rest of the spins located on the same 
side with respect to the impurity. In summary, what is 
predicted is a very local enhancement of the correlations 
around the spin-1 impurity. The region involved in the 
process corresponds to only 3 spins (the spin-1 and the 
two neighboring spin- 1/2). From the point of view of 
the rest of the spins on the chain, this three-spin struc- 
ture has formed a singlet and it behaves effectively as 
a nonmagnetic vacancy. Thus, this result for dimerized 
chains agrees with a general result obtained for uniform 
AF Heisenberg chains which states that an impurity with 
spin-S is equivalent to open boundary conditions with a 
decoupled spin of size (S-l). JnJ Hence, a spin-1 impu- 
rity is equivalent to a non-magnetic impurity and forms 
a singlet. 

These features can be clearly seen in Table |. The 
correlations S(l) have been computed with a standard 
world-line Monte Carlo at a temperature T = 0.0625J 
and for TV — 80, except otherwise stated. The dimeriza- 
tion is S = 0.05. The statistical error of the QMC simula- 
tions affects only the last digit of the quantities reported. 
The effects of nonmagnetic impurities are shown for the 
40 site chains for the case where OBC are used. Results 
for the 39 site chains are very similar to the values for the 
40 site chain. The effects of 5 = 1 impurities are shown 
in the last four rows. In these cases the position of the 
second impurity iq is shown (the first one is located at 
site 1). In all cases, it has been indicated if the bond 
between the reference site, i\, and its spin-1/2 nearest 
neighbor site is strong (s) or weak (w). In Table [ijj, the 
dependence of the previous results with 5 is studied using 
Lanczos diagonalization results for a chain with TV = 20 
sites and one impurity. Note the fairly good agreement 
between these zero temperature results compared with 
the finite temperature QMC results on the 80 site chain 
both for the pure system and for the case of two impu- 
rities located at sites io = 1,41 and 5 — 0.05. It can 
be observed on Table || that as S is increased from 0.05 
to 0.5 there is a smooth variation in the behavior of the 
various correlations shown. In all cases 5(1), 5(2) and 
5(3) are reduced from their pure values, and 5j mj) (l) and 

Across 

(2) are enhanced. 
The possibility raised above that the singlet ground 
state of the three-site subsystem 1/2-1-1/2 could be- 
have as an effective nonmagnetic vacancy, and hence 
that it may lead to an enhancement of the AF corre- 
lations |2l],|22| between the sites close to this vacancy, is 
actually confirmed by the results shown in Tabic III. In 
this Table, the correlations between a spin next to one of 
these three-site clusters, i.e. located at two lattice spac- 
ings from the 5=1 impurity, and the remaining spins of 
the chain are shown to be enhanced with respect to their 
values in the pure system. As previously noticed |^l|,^2| , 
the enhancement of the correlations for the case when 
the bond between the two first sites next to a vacancy is 
weak is more pronounced than the case of a strong bond 



between them. 

TABLE I. Spin-spin correlations S(l) at various distances 
as denned in the text obtained by Monte Carlo simulations on 
chains with TV = 80 sites, using periodic boundary conditions, 
and at T = 0.0625 unless otherwise stated. "Pure" means 
that there are no impurities. The chains with TV = 39, 40 sites 
have open boundary conditions, and the correlations are mea- 
sured starting at the site next to the boundary. The dimer- 
ization is 8 — 0.05. The last four rows correspond to the chain 
in the presence of two 5 = 1 impurities, the first one located 
at site 1 and the second one at site io- The reference site is lo- 
cated at ii. w, s refer to weak and strong bonds, respectively, 
between the spin at site ii and its spin-1/2 neighbor. 





S(l) 


S(2) 


S(3) 


Across (2) 


Simp (l) 


pure s 


-0.777 


0.230 


-0.247 






pure s (T=0.05) 


-0.785 


0.230 


-0.250 






pure w 


-0.400 


0.230 


-0.122 






pure w (T=0.05) 


-0.408 


0.230 


-0.126 






N=40 s 


-0.906 


0.242 


-0.274 






N=40 s (T=0.05) 


-0.907 


0.244 


-0.275 






N=40 w 


-0.777 


0.297 


-0.253 






N=40 w (T=0.05) 


-0.777 


0.299 


-0.254 






imp. io=40 ii=41 s 


-0.408 


0.170 


-0.097 


0.295 


-0.971 


imp. io=40 ii =39 w 


-0.313 


0.189 


-0.081 


0.295 


-1.067 


imp. io=41 ii=40 s 


-0.429 


0.196 


-0.141 


0.293 


-1.148 


imp. io=41 ii =42 w 


-0.293 


0.193 


-0.087 


0.293 


-1.240 



TABLE II. Spin-spin correlations at various distances as 
defined in the text obtained by Lanczos diagonalization for 
TV = 20, using periodic boundary conditions and for several 
values of the dimerization <5. 





S(l) 


S(2) 


S(3) 


Across (2) 


Simp ( 1) 


pure w (8 = 0.05) 


-0.4113 


0.2319 


-0.1339 






pure w (8 = 0.1) 


-0.3213 


0.2111 


-0.0927 






pure w (8 = 0.3) 


-0.1706 


0.1401 


-0.0329 






pure w (8 = 0.5) 


-0.0954 


0.0860 


-0.0117 






pure s (8 = 0.05) 


-0.7679 


0.2319 


-0.2466 






pure s (8 — 0.1) 


-0.8459 


0.2111 


-0.2378 






pure s (8 = 0.3) 


-0.9503 


0.1401 


-0.1574 






pure s (8 = 0.5) 


-0.9834 


0.0860 


-0.0931 






w (8 = 0.05) 


-0.3012 


0.1953 


-0.0976 


0.2936 


-1.2385 


w (8 = 0.1) 


-0.2579 


0.1828 


-0.0760 


0.2838 


-1.2600 


w (8 = 0.3) 


-0.1511 


0.1272 


-0.0300 


0.2179 


-1.3018 


w (8 = 0.5) 


-0.0870 


0.0795 


-0.0110 


0.1408 


-1.3211 


s {8 = 0.05) 


-0.4325 


0.1951 


-0.1438 


0.2936 


-1.1462 


s (8 = 0.1) 


-0.5135 


0.1832 


-0.1569 


0.2838 


-1.0723 


s (8 = 0.3) 


-0.7789 


0.1315 


-0.1396 


0.2179 


-0.7346 


s {8 = 0.5) 


-0.9206 


0.0841 


-0.0901 


0.1408 


-0.4342 



G 



TABLE III. Spin-spin 

correlations 5(2, Z) = (5 Z (2)5 2 (2 + Z)) at various distances 
from a site two lattice spacings apart from a 5 = 1 impurity 
obtained by Lanczos diagonalization for N — 20 for various 
values of 8, and by Monte Carlo simulations for iV = 80 sites 
and 8 = 0.05. 





S(2,3) 


S(2,4) 


S(2,5) 


S(2,6) 


pure s (5 = 0.05) 


-0.7679 


0.2319 


-0.2466 


0.1243 


pure s (8 = 0.10) 


-0.8459 


0.2111 


-0.2378 


0.0960 


pure s (8 = 0.30) 


-0.9503 


0.1401 


-0.1574 


0.0360 


pure s (5 = 0.50) 


-0.9834 


0.0860 


-0.0931 


0.0125 


s (S = 0.05) 


-0.8224 


0.2371 


-0.2580 


0.1255 


s (8 = 0.10) 


-0.8696 


0.2141 


-0.2421 


0.0976 


s (5 = 0.30) 


-0.9531 


0.1404 


-0.1577 


0.0361 


s (S = 0.50) 


-0.9838 


0.0861 


-0.0931 


0.0125 


pure w (5 = 0.05) 


-0.4113 


0.2319 


-0.1339 


0.1243 


pure w (5 = 0.10) 


-0.3213 


0.2111 


-0.0927 


0.0960 


pure w (<5 = 0.30) 


-0.1706 


0.1401 


-0.0329 


0.0360 


pure w (8 = 0.50) 


-0.0954 


0.0860 


-0.0117 


0.0125 


w (5 = 0.05) 


-0.6263 


0.2708 


-0.2113 


0.1602 


w (8 = 0.10) 


-0.5069 


0.2633 


-0.1606 


0.1417 


w (5 = 0.30) 


-0.2183 


0.1678 


-0.0464 


0.0492 


w (5 = 0.50) 


-0.1040 


0.0926 


-0.0133 


0.0141 


N=80 pure s 


-0.777 


0.230 


-0.247 


0.115 


N=80 pure w 


-0.400 


0.230 


-0.122 


0.116 


N=80 i =40 s 


-0.815 


0.231 


-0.251 


0.117 


N=80 i =40 w 


-0.619 


0.260 


-0.200 


0.143 


N=80 i =41 s 


-0.829 


0.235 


-0.260 


0.120 


N=80 i =41 w 


-0.629 


0.271 


-0.210 


0.157 




FIG. 3. Dynamical structure factor S zz (q,uj) obtained by 
exact diagonalization in the 16 site chain with two 5=1 im- 
purities for 8 = 0.05 averaged over all the impurity positions. 
The energy of the first peak for the pure N = 16 chain are 
indicated with full circles. 



In summary, the results shown thus far seem to in- 
dicate that for magnetic impurities the AF correlations 
between these impurities and their two neighboring spin- 
1/2's are strong. However, note that once the three-spin 
cluster centered at the impurity is considered as a singlet, 
then the rest of the spins will enhance the correlations 
among themselves as if in the presence of a nonmagnetic 
impurity. 

Finally, in Fig.| S zz (k ,u>) is shown for N = 16, 
5 = 0.05 and in the presence of two 5 = 1 impurities. 
The results are averaged over all the impurity positions. 
The energies of the first peaks for the pure system have 
been indicated for reference. The behavior of S zz (k,u>) 
in the presence of S = 1 impurities is very similar to 
the previously reported behavior for nonmagnetic impu- 
rities p5| . In other words, two bands appear here: 1. 
the low energy "impurity band" with the highest inten- 
sity peaks at or near k — tt, and 2. a high energy band, 
which is essentially a remnant of the dispersion relation 
for the undoped system. 



IV. MAGNETIC SUSCEPTIBILITY AT FINITE 
TEMPERATURES 

The suppression of the spin gap due to magnetic impu- 
rities can also be detected in thermodynamical properties 
such as the magnetic susceptibility x(T). It is well-known 
that the presence of a spin gap between the singlet ground 
state and the excited triplet state leads to an exponential 
decay as the temperature is reduced to zero. On the other 
hand, x(T) would diverge or converge to a nonzero value 
as T vanishes for a gapless system, such as the uniform 
ID Heisenberg model. A similar result should emerge as 
a consequence of the presence of states inside the original 
gap induced by impurities, as discussed in the previous 
section. This change of the magnetic susceptibility upon 
Zn, Mg or Ni doping was indeed observed experimentally 

ijbj^flG 44 In particular, in 



in the SP compound CuGeOs 
susceptibility measurements along the c axis it has been 
reported the presence of a second maximum at temper- 
atures lower than the SP transition temperature. |L4 44 
This second maximum in the susceptibility along the c 
axis has been attributed to the transition to a long-range 
AF order and it occurs at a temperature slightly higher 
than the corresponding Neel temperature (Tn)- 

In Fig.^, the magnetic susceptibility computed by 
QMC simulations is shown as a function of temperature 
for a dimerized Heisenberg chain with 8 = 0.05 in the 
presence of both magnetic and non-magnetic impurities. 
The statistical error of the QMC simulations is about the 
size of the open circles in this figure. The susceptibility 
for the pure system on a chain of 80 sites has also been 
computed. The results for this case are very similar to 
the case of N = 40 with OBC and strong edge bonds and 
they are not included in the figure for clarity. %(T) for 
N = 80 with two spin-1 impurities has also been com- 
puted, for the case where the first impurity is located at 
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site 1 and the second one at various sites such as 20, 21, 
40 and 41. The susceptibilities obtained for these systems 
are compared with those obtained from chains of N — 40 
with open boundary conditions. The latter chain can be 
considered as resulting from a N — 80 chain after a pair 
of nonmagnetic impurities have been introduced. 

In the case of N = 40, OBC, two possibilities have been 
also analyzed: the case where the bonds at the chain ends 
are both weak and the case where they are both strong. 
The results can be interpreted in terms of the intuitive 
picture developed in Ref. | p5[ for the large dimerization 
limit. In the case of strong end links, the spins are basi- 
cally paired in strong dimers and hence it is quite likely 
that a spin gap is still present in the chain leading to 
an exponential decay of x(T) as T — > observed in the 
figure. This behavior is not observed in the other case of 
weak end bonds where there are loose spins at the ends 
of the chain which cause a strong suppression of the spin 
gap. In the case of spin-1 impurities, depending on their 
relative positions, a periodic chain is effectively cut, or in 
chains with odd number of sites, or in chains with even 
number of sites where at least one of them has both weak 
end bonds. In both cases, the susceptibility would not 
decay exponentially to zero as the temperature goes to 
zero as it is observed in the figure at least down to the 
lowest temperature achievable with our algorithm. 




0.00 



0.0 



0.2 



0.4 



0.6 



T/J 



FIG. 4. Magnetic susceptibility as a function of temper- 
ature obtained by Monte Carlo simulations in dimerized 
Heisenberg chains with S — 0.05. N — 80, two S = 1 im- 
purities in sites 1 and 40 (up triangles), in sites 1 and 41 
(open circles), in sites 1 and 20 (filled diamonds), and in sites 
1 and 21 (open diamonds). TV = 40, OBC, weak end bonds 
(stars), strong end bonds (crosses). In the inset, it is shown 
X(T) obtained by exact diagonalization in chains of N = 12 
(solid line), TV = 14 (short dashed line), and N = 16 (long 
dashed line), and OBC with weak end bonds. The results ob- 
tained by QMC for TV = 40, OBC, weak end bonds (circles) 
for temperatures down to T/J = 0.04 are also shown. 



In principle, the experimentally observed behavior of 
the susceptibility should be theoretically accounted for by 
averaging the results obtained for chains of even number 
of sites with the ones obtained for chains with odd num- 
ber of sites. The latter diverge as T — > due to the total 
S z = 1/2. For the case of nonmagnetic impurities (open 
boundary conditions), from both QMC and exact diago- 
nalization results for chains up to 16 sites, we conclude 
that by averaging the susceptibilities of chains with even 
number of sites and strong end bonds with the suscepti- 
bilities of chains with odd number of sites, one does not 
obtain a second maximum in the low temperature region. 
The Curie behavior of the odd chains dominates over the 
exponentially decaying behavior of the even chains at low 
temperatures. On the other hand, the susceptibility on 
chains with even number of sites and weak end bonds, 
as it can be seen in the inset of Fig.^, does indeed show 
this second maximum and it is quite likely that this sec- 
ond maximum can survive the averaging procedure as we 
have verified numerically. 

Let us examine to what extent these results are rele- 
vant to spin-Peierls systems. In these systems, upon the 
introduction of nonmagnetic impurities, the spin-lattice 
coupling leads to the case of chains where the two edge 
bonds are strong, even for an odd total number of sites in 
which case a soliton is formed in the middle of the chain. 
PU However, it has been recently claimed that elastic 
interchain interactions may stabilize a given pattern of 
distortions, i.e. with and without nonmagnetic impuri- 
ties the weak and strong bonds of the dimerized chain 
may remain the same as for non-dynamical phonons flsj ] 
and this behavior has been verified numerically. |3i| That 
is, even site chains with both weak links are still present 
in SP chains due to elastic coupling to pure neighbor 
chains. 

Besides, it is reasonable to think that the results ob- 
tained for larger open chains are going to be indicative of 
smaller nonmagnetic impurity doping. Then, as it can be 
seen in the inset of Fig.|| the shift to lower temperatures 
of the second peak as N is increased is in complete agree- 
ment with the experimental results showing that the Need 
temperature decreases as the impurity concentration de- 
creases 0] . One should notice that this behavior of Tn 
occurs in the dimerized region of CuGeOs where precisely 
our fixed-dimerization study should apply. (A different 
behavior in the uniform region for large impurity concen- 
trations has been reported H]). Finally, the fact that the 
second maximum in x is due to the contribution of even 
site chains with weak end links is very consistent with 
the fact that, as mentioned above , these chains present 
the largest enhancement of AF correlations upon doping. 

On the other hand, the susceptibility in the case of 
two magnetic impurities has already a very similar tem- 
perature dependence for different impurity positions and 
it is quite likely that an average over all these impurity 
positions would still result in a curve of the same shape. 
Based on the mapping of this problem to an effective 
nonmagnetic impurity chain a second maximum should 
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exist in \ an d it should appear at temperatures lower 
than the attainable ones within our QMC algorithm. 



V. CONCLUSIONS 

In this article the effects of S = 1 magnetic impuri- 
ties on dimerized AF Heisenberg chains have been an- 
alyzed, and the results have been compared with those 
originated by the introduction of nonmagnetic impurities. 
The overall effects, such as the destruction of the spin gap 
which reveals itself in the dynamical spin structure factor 
S(q, uj) (Figs. || and||) and in the magnetic susceptibility 
(Fig. 0), are similar between the two types of impurities. 
However, the microscopic origin and details of the tight 
states formed between the impurities and the spin-1/2 
sites close to it are different. This difference translates 
into the existence of strong AF correlations in the imme- 
diate vicinity of the S = 1 impurity. This local behavior 
could eventually be experimentally detected, for example 
by using NMR techniques. On the other hand, the spins 
located at distances of two lattice spacings (or larger) 
from the impurity react to the spin singlet state aris- 
ing from the coupling of the S=l impurity with its two 
neighboring spin-1/2 as if it were a nonmagnetic impu- 
rity. This interesting effect produces an enhancement of 
AF correlations in the vicinity of the S = 1 impurities as 
for Zn doping, which may lead to a global stabilization of 
the AF state as found experimentally in CuGe03 using 
several techniques. 

As discussed in the previous section with respect to 
the case of nonmagnetic impurities, it is expected that in 
spin-Peierls systems the partial substitution of Cu ions 
by Ni ions would also change the dimerization pattern 
which these compounds present below the SP transition 
temperature. That is, the changes in the spin-spin cor- 
relations discussed in Section HI would imply, due to the 
spin-phonon coupling, a distortion of the dimerization 
pattern, especially near the Ni impurities. Hence our 



results could in principle not be directly related to the 
experimentally observed behavior of Ni-doped CuGe03. 
However, some previous studies, as well as new results 
which we have recently obtained, indicate that the elas- 
tic coupling between a doped chain with its neighbor un- 
doped ones tends to stabilize the undistorted dimerized 
pattern and in this case our results could be directly rel- 
evant to experiments in SP systems. In this sense, the 
agreement of our results for the magnetic susceptibility 
with recent experimental results on CuGeOa is remark- 
able and this leads to an interesting explanation on the 
origin of the second peak observed in these experiments. 
On the other hand, the effects of S = 1 impurities ana- 
lyzed in this paper could be also detected in the struc- 
turally dimerized compounds such as the ones indicated 
in the Introduction. 
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